The (3^ potential at the U(5)— 0(6) critical point of the Interacting Boson Model 
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Exact numerical results of the interacting boson model Hamiltonian along the integrable line 
from U(5) to 0(6) are obtained by diagonalization within boson seniority subspaces. The matrix 
Hamiltonian reduces to a block tridiagonal form which can be diagonalized for large number of 
bosons. We present results for the low energy spectrum and the transition probabilities for systems 
up to 10000 bosons, which confirm that at the critical point the system is equally well described by 
the Bohr Hamiltonian with a f3'^ potential. 
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I The main goal of this note is to report on new results that complete a previous study JJ on the relations between 
, ' the critical point in the transition from U(5) to 0(6) limits of the interacting boson model (IBM) and the recently 
^ ' proposed E(5) critical point symmetry Q. In Ref. |[1| two different boson Hamiltonians performing the transition 
from U(5) to 0(6) were used to show that at the critical point: i) they provide different spectra and transitions for 
' small number of bosons, ii) they converge to the same spectrum for large N and iii) both converge to the spectrum 
, provided by the Bohr Hamiltonian with a /J^ potential rather than to the one provided by a square well potential 
as in the E(5) model. Most of the large N analysis was based on the solution of the Richardson equations 
which allow to obtain energy eigenvalues but the form of the eigenstates is not well suited to calculate transition 
probabilities. Therefore, in order to study transition rates we had to resort to current IBM codes that restricted 
our calculations to systems up to N=40 1]. For these small N values, the transition rates show a tendency to 
approach the (3^ potential results but they are not conclusive. In this brief report we present an alternative to the 
Richardson equations for obtaining the low energy eigenvalues and, on the same footing, the transition probabilities 
. in the U(5)-0(6) transitional region for large N values. 

' The two Hamiltonians describing the U(5)-0(6) transition studied in Ref.m are 
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Hi = xhd + j;^P^P , (1) 
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Hii = xfid - ^0^=° • 0^^=° , (2) 



where 

fid 



Q'^=° = (fit xd + d^ X S)(2), (5) 

and • stands for the scalar product. We have introduced in Q the boson pair creation operators P| = • and 
= gt . i^jiat will be used later on. 
The mean field analysis of the quantum phase diagram of the IBM is usually performed within the intrinsic state 
formalism 0)13 where, after separating the three Euler angles, the trial wave function is a boson condensate depending 
on the two geometrical variables f3 and 7. Along the U(5)-0(6) transition the energy surface is 7— independent and the 
intrinsic ground state energy for a given value of the control parameter x corresponds to the value of the deformation 
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parameter, /?, which minimizes the energy surface. The phase transition along this Une is then signaled by the 
condition 

[d^E{N,P)/dp\=o^Q, (6) 

which fixes the critical value of the control parameter x. For the Hamiltonian the critical x\sx^^ — 0.5, independent 
of the number of bosons iV, while for the Hamiltonian (O it is x^J — |^5f • In the large N limit x[^ — > 4/5. 

In order to study the eigenstates of the Hamiltonians ^ and Q we introduce the s— and d— boson pair algebra 
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For each £ value, or 2, the three operators {K'^ ^ , K^} satisfy the su(l, 1) commutator algebra 

[K^,K^] = ±6ee,Kf^ , [K+,K-] = -2Su'K^. (9) 

A complete set of eigenstates for a general IBM U(5)-0(6) transitional Hamiltonian can be written in terms of the 
raising operator K'^ acting on a subspace of unpaired bosons characterized by the seniority quantum number ly 
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\nm)^^={K+Y'\vt), (10) 



where = 0, 1, i/^ = 0, 1, 2, . . . , and \vi) is a normalized state. The value of vg gives the number of bosons of type I 
not coupled in pairs to zero angular momentum. The label hi refers to boson pairs coupled to zero angular momentum. 
Therefore, the total number of bosons is iV = 2hs + 2hd + + Vd- Using the su(l, 1) algebra it is straightforward to 
obtain the normalization constants 

r^"* I /'^r-rVs'+rM,, \ ni\{2ht, + 2i^2vi-\)\\ 

Now we proceed to construct the complete set of states as 

\nshd, VsVd) = = (A'+)"° (-ftrj")"" WsVd) ■ (11) 



The basis (|ll|l . although it has not good angular momentum, is especially useful for diagonalizing the Hamiltonians 
and ISJ. To show this, we rewrite the Hamiltonians Q and ((5J in terms of the generators of the two su(l, 1) 
algebras 

iii = + ^1-^ {K^ + X+X,- - Kt Kd - K+R-) , (12) 

Hij = xhd-^-^{4KtK^ +iK+K; +5h, + 7id + 2hshd) . (13) 
The matrix elements of the relevant operators for both Hamiltonians in the basis (|ll|l are: 

{hshd,VsVd\hs\hshd,VsVd) = 2hs + Vs 
{hshd,VsVd\nd\nshd,VsVd) = 2hd + Vd 

{hsfid, I's'^dl \hshd, i^si^d) = hdihd + Vd + 
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{{hs-l){hd + l) ,VsVd\KtK^\hshd,VsVd) = -^^s (fid + 1) (2ns + '2vs - 1) (2^^ + 2i^d + 5). (14) 
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The Hamiltonians and jSJ do not mix states with different seniority quantum numbers {vg, Vd) leaving invariant 
these seniority subspaces. Within each subspace the Hamiltonian matrices are tridiagonal and can be easily diago- 
nalized for very large N values. We will label states within each subspace by the quantum number ^. It is worthwhile 
to note here that d boson seniority, is equivalent to the 0(5) quantum number t 2]. The construction of the 
spectrum for a system with even number of bosons is as follows, one starts with the subspace t = [vg — Q^Vd — 0) 
where all the bosons are coupled in pairs of zero angular momentum. Consequently, states within this subspace will 
have total angular momentum L — Q. The lowest eigenvalue (^ = 1) is the ground state O^q (the notation is L^^ 
0), the second lowest eigenvalue is the first excited state r = L'^ = 0"^ state, which is labelled OJq, etc. The next 
block with T — \ {vs — \,Vd — 1) has one pair broken into an s-boson and a c?-boson. Correspondingly, all states 
in this block have L = 2. The lowest eigenvalue (^ = 1) is the lowest 2+ which is labelled as 2^^^ the next one 
(^ = 2) is 22 1, 6tc. The next block is for r = 2 [ug — Q^Vd = 2) and corresponds to one broken boson pair into two 
d-bosons. It provides states with angular momenta L = 4, 2 (notice that L = is excluded from this subspace since it 
is included in the = 0, z^d = subspace). One can continue in this way with the next block, t = Z [vg ^ l,Vd = 2>). 
It corresponds to two boson pairs broken into one s— boson and three d— bosons and gives rise to L = 6, 4, 3, states 
and so on. The ground state band is formed by all lowest (^ — 1) eigenstates for r — 0, 1, 2, 3, . . . The first excited 
band (i^ — 2) is formed by the next lowest t = 0, 1, 2, 3, . . . eigenstates, etc. Following this sequence one finds the well 
known triangular structure associated to 0(5). All this is shown schematically in Fig.^ 
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FIG. 1: Schematic spectrum obtained by diagonalization within boson seniority subspaces as explained in the text and its 
correspondence with the one of Refs. 0] and 

The diagonalization of the Hamiltonian in each subspace provides the necessary information to calculate electro- 
magnetic transition rates. We will be interested here on the electric quadrupole transitions, which apart from an 
unimportant scale factor, are described by the quadrupole operator ij^l. The action of this operator on the basis 
states without broken pairs is 



Q/.|^is"-d,0, 0) 



^s^iS^d.Q 



{Ktr-' {Kjf^ 10) + ud {Ktr (Kir^' m 



(15) 



where |0) stands for the boson vacuum. 

The matrix elements of interest if one wants to evaluate transition rates from the ground state to the first excited 
state are 



{{hg - 1) hd, 1, 1| \nghd, 0, 0) 



2fig (2nd + 5) 



(16) 
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(hs [fid - 1) , 1, 1| Q^i |nsn.d,0, 0) = 



2nd (2ns + 1) 



(17) 



If we write the eigenstates as 



I*, VsVd) = ^ \nsnd, VsVd) , 



(18) 



the matrix element of the Q operator between the gromid state, Y^, 00), and the first excited state, |$, 11), is 

4 
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FIG. 2; Variation with the number of bosons (up to = 40) of selected energy and B(E2) ratios for IBM calculations performed 
at the critical points of Hamiltonians Q (broken line) and ^ (full line). The corresponding E(5) and values are marked 
with horizontal dotted lines. 



($,ll|QJvi/,00)= 



2n, {2nd + 5) 



.00 



2hd {2hs + 1) 00 11 



(19) 



The matrix elements of the electric quadrupole operator between the first excited state {vg = \,Vd = 1) and the states 
with Us = Q,Vd = 2 can be calculated in a similar way. 

In Fig. 2 we present some selected low energy eigenvalues and B{E2) ratios for boson numbers up to — 40 at 
the critical points of both IBM Hamiltonians JQ) and We would like to emphasize here that the critical points 
for the two Hamiltonians are different. The four energy ratios presented are written explicitly in the figure and the 

two displayed B(E2) ratios are: Ri — ' Y — 7pr^ and i?2 — , ' — 4^. Where we are using the notation 

to indicate the states. The purpose of this figure is to correct a mistake we had in Fig. 3 of Ref. |l| where the results 
for the Hamiltonian (01 were calculated with a wrong value for Xc- As can be seen in the figure, there are sizable 
differences in the spectrum and transition rates between both Hamiltonians at the critical points. Though from Fig. 
2 a general tendency for convergence to the solution of the Bohr equation with a Z?** potential rather than to the E(5) 
symmetry is observed, the results, especially from the B(E2)'s, are not yet conclusive. In Fig. 3 we show the new 
results of this report. The same quantities as in Fig. 2 are plotted for N values up to 10000, including transition 
rates. In Ref. tl| energy eigenvalues were calculated up to = 1000 and transition rates up to = 40. We can 
now clearly appreciate, both from the energies and B(E2) transitions, that the IBM Hamiltonians at the U(5) to 0(6) 
critical point in the large N limit converge to the Bohr Hamiltonian with a /3'* potential. 
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FIG. 3: Same as Fig. |5|but here the number of bosons runs up to 10000 in both the energy and the B(E2) ratios. Please note 
the logarithmic scale in the number of boson axis. 



In this report we make use of the property that the IBM Hamiltonian along the transitional line from U(5) to 0(6) 
is block diagonal wit respect to the boson seniority quantum number and tridiagonal within each subspace. This 
reduction allows to obtain exact solutions up to very large number of bosons for energies and wave functions. We 
have applied this formalism to confirm previous studies about the correspondence between the IBM Hamiltonians 
at the critical point in the U(5)-0(6) transition and the Bohr Hamiltonian with a potential for the low energy 
properties. This issue has also been studied recently from a different point of view by Rowe et al. [To|. The formalism 
presented here for the IBM can be easily generalized to other two- level boson models 11]. 
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